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ABSTRACT 

Aims. We study the proximity effect in the Lyo- forest around high redshift quasars as a function of redshift and environment employing 
a set of 3D continuum radiative transfer simulations. 

Methods. The analysis is based on dark matter only simulations at redshifts 3, 4, and 4.9 and, adopting an effective equation of 
state for the baryonic matter, we infer the HI densities and temperatures in the cosmological box. The UV background (UVB) and 
additional QSO radiation with Lyman limit flux of L,, LL = 10 31 and 10 32 erg Hz" s" 1 are implemented with a Monte-Carlo continuum 
radiative transfer code until an equilibrium configuration is reached. We analyse 500 mock spectra originating at the QSO in the most 
massive halo, in a random filament and in a void. The proximity effect is studied using flux transmission statistics, in particular with 
the normalised optical depth £ = T e ff Qso/T"cff,Lytt, which is the ratio of the effective optical depth in the spectrum near the quasar and 
in the average Lyo- forest. 

Results. Beyond a radius of r > 1 Mpc h from the quasar, we measure a transmission profile consistent with geometric dilution 
of the QSO ionising radiation. A departure from geometric dilution is only seen, when strong absorbers intervene the line-of-sight. 
The cosmic density distribution around the QSO causes a large scatter in the normalised optical depth. The scatter decreases with 
increasing redshift and increasing QSO luminosity. The mean proximity effect identified in the average signal over 500 lines of 
sight provides an average signal that is biased by random large scale density enhancements on scales up to r » 15 Mpc h~ l . The 
distribution of the proximity effect strength, a parameter which describes a shift of the transmission profile with respect to a fiducial 
profile, provides a measure of the proximity effect along individual lines of sight. It shows a clear maximum almost without an 
environmental bias. Therefore this maximum can be used for an unbiased estimate of the UVB from the proximity effect. Differing 
spectral energy distributions between the QSO and the UVB modify the profile which can be reasonably well corrected analytically. 
A few Lyman limit systems have been identified that prevent the detection of the proximity effect due to shadowing. 
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1. Introduction 
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The baryon content of intergalactic space is responsible for the 
large number of absorption lines observed in the spectra of high 
redsh ift quasars. This phenomenon, al so known as the Lya for- 
est dSargent etalj [l980; Rauch 1998), is mainly attributed to 
intervening Hi clouds along the line of sight (LOS) towards a 
QSO. The majority of the Hi Lya systems are optically thin to 
ionising radiation. The gas is kept in a high ionisation state by 
the integrated emission of ultra-violet photons originating from 
the overall populations of quasars an d star-forming galaxies : 
the ultra-violet ba c kground field (UVB: lHaardt & Madaulll996t 
iFardal et all 1 19981: lHaardt & Madaull200ll) . Accurate estimates 
of the UVB intensity at the Lyman limit are crucial for the un- 
derstanding of the relative contribution of stars and quasars to 
the UVB and also to ensure realistic inputs for numerical simula - 
tions of structure formation dDave et alJl999HHoeft et al.l2 006). 

The intensity of the UVB can be sttongly affected by en- 
ergetic sources such as bright QSOs leading to an enhance- 
ment of UV photons in its environment. The neutral fraction 
of H i consequently drops yielding a significant lack of absorp- 
tion around the QSOs within a few Mpc. Knowing the lumi- 
nosity of the QSO at the Lyman limit, this "proximity effect" 
has been widely used in estimating the i ntensity of the UVB 
dCarswell et all Tl98 7: Baitli k et alfll988l) at various redshifts 
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on large sample s (Baitli k et alll988l:lLu et alj|1991t [Scott et al.l 
2000; Liske & Willigerl 1200 lb IScott et alJ 120021 but compare 
iKirkman & Tvtlerj 120081) and recently also towards individual 
lines of sight dLu et alJI 19961: iDall'Aglio et alj|2008bl) . 

Principal obstacles in the analysis of the proximity effect ar- 
rise due to several poorly understood effects which in the end 
might lead to biased estimates of the UVB. The most debated 
influence is the gravitational clustering of matter around QSOs. 
Clustering enhances the number of absorption lines on scales of 
a few proper Mpc, which if not taken into account might lead to 
an overestimation of the UVB by a factor of a few. Alternatively, 
star forming galaxies in the QSO-environment may contribute 
to an enhancement, or over-ionisation, of the local ionised hy- 
drogen fraction, i.e. we would underestimate the UVB from the 
measured QSO luminosity. In the origin al formalism o f the prox- 
imity effect theory introduced by Baitl ik et alJ JT988) (hereafter 
BDO), the possibility of density enhancements near the QSO 
redshift was neglected . Howe ver, lLoeb & Eisensteinl (1 19951) and 
iFaucher-Giguere et al.l ([2008) found that the UVB may be over- 
estimated by a factor of about 3 due to high en vironmental den- 
sity. C omparing observations and simulations, iD'Odorico et al.l 
(2008) found large scale over densities over about 4 proper Mpc. 
The disagreement between the UVB obtained via the proximity 
effect and from flux transmission statistics has be en used to esti- 
mate the average density profile around the QSO dRoll inde et al. 
2005; Guimaraes et al 1 l2007tlKim & Croftll2008l) . Only recently 
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iDall'Aglio etai](l2008ah showed that a major reason for overes- 
timating the proximity effect is not a general cosmic overdensity 
around QSOs, but the methodological approach of estimating 
the UVB intensity by combining the proximity effect signal over 
several sight lines. Providing a definition of the proximity effect 
strength for a given LOS, they propose instead to investigate the 
strength distribution for the QSO sample, yielding consistency 
with the theoretical estimates of the UVB. 

Numerical simulations of structure formation have been 
a crucial tool in understanding the nature and evolution of 
the Lya forest. One of the major challenges in the current 
development of 3D simulations represents the implementa- 
tion of radiative tr ansfer into the forma l ism of hierarchical 
structure formation dGnedin & Abel 1200 lb iMaselli et all [2 003 ; 



Razoumov & Cardall 2005; Riikhorst et al. 2006; Mellema et al. 
2006; iPawlik & Schavd 120081) . In particular the coupling of 
radiative transfer with the density evolution requires a large 
amount of computational resources and remains up to now a 
largely unexplored field. Therefore, most available codes apply 
radiative transfer as a post processing step. 

Only a couple of radiative tra n sfer studies exist for th e 
Lya forest dNakamoto et all 1200 it iMaseUi & Ferraral [2005). 
They discuss the influence of radiative trans fer on the widel y 
used semi-analytical mo del of the forest by iHui et al] d!997l) . 
IMaselli & Ferra ra (2005) find that due to self-shielding of the 
UVB flux in over-dense regions the hydrogen photo-ionisation 
rate fluctuates by up to 20 per cent and the helium rate by up 
to 60 per cent. Since radiative transfer is not negligible for a 
background field, this should be true for point sources as well. 
Up to now numerical studies dealt with ionisation bub bles in the 
pre-reionisation era or right at the end of reionisation ([iliey et al] 



pre-reionisation era or right at the end of reionisation (lliev et al. 
2007tlKohleret al] l2007: Ma selli et aljl2007HTrac & CerJl2007t 



Zahn et al. l2007ir These studies aim at determining the sizes of 



H ii regio ns as observed in transm ission gaps of Gunn-Peterson 
troughs dGunn & Peterson!! 19651) in high-redshift sources, or in 
21cm emission or absorption signals. By comparing radiative 
transfer simulations with synthetic spectra, Mas elli et al] (12007!) 
showed, that the observationally deduced size of the H n regions 
from the Gunn-Peterson trough is underestimated by up to 30 
per cent. 

In the study of the proximity effect at redshifts lower than 
that of reionisation, the influence of radiative transfer has been 
neglected up to now, as the universe is optically thin to ionising 
radiation. We expect that absorption of QSO photons by inter- 
vening dense regions might reduce the QSO flux. These dense 
regions can shield them selves from the QSO radiation field 
dMaselli & Ferr ara 2005) and would not experience the same 
increase in the ionisation fraction as low density regions. This 
leads to a dependence of the proximity effect on the QSO envi- 
ronment. Furthermore, the amount of hard photons in the QSO 
spectral energ y distribution affects the p roximity effect profile as 
suggested bv iDalF Aglio et al . (2008b). Similar to the study of 
high redshift H n regions where ionisation fronts are broadened 
due t o the ionising flux' s shape of the s pectral energy distribu- 
tion dShapiro et aljl2004lOiu et al.ll2007h . we expect the size of 
the proximity effect zone to be a function of the spectral hard- 
ness. Harder UV photons have smaller ionisation cross sections 
and cannot ionise hydrogen as effectively as softer UV photons. 

In this study we employ a three dimensional radiative trans- 
fer simulation to study H n regions expanding in a pre-ionised 
intergalactic medium (IGM) at redshifts z = 3,4, and 4.9. We 
intend to quantify the influence of the above mentioned ef- 
fects. To this end, we consider two realistic QSO luminosities. 
Furthermore, we study the influence of environment by placing 
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Fig. 2. Mock spectra synthesised from the z — 3 snapshot with a 
QSO in a filament. The QSO sits on the right hand side of each 
spectrum, and the wavelength is in Angstrom. The dashed line 
gives the Lya forest spectrum without the influence of the QSO, 
the solid line the spectrum from the radiative transfer simulation 
of the proximity effect for a QSO with a Lyman limit luminos- 
ity of L VhL = 10 31 erg Hz s _1 . The spectra have a resolution of 
Av = 6.7 km s~ and a signal to noise of 70. 



the QSO either inside a massive halo, in a random filament, and 
in an under-dense region (void). 

The paper is structured as follows. In Sect.[2]we discuss the 
dark matter simulation used in this study, and we show that a 
realistic model of the Lya forest is found from a semi-analytical 
model of the IGM. In Sect. [3] we describe the method to solve 
the radiative transfer equation, and we discuss in Sect. 13.21 and 
!3.3l how the different UV sources are implemented. In Sect. [4] we 
review the standard appro ach used to characterise the proximity 
effect dBaitlik et al.ll 1988b . In Sect. [5] we describe the different 
radiative transfer effects on the overionisation profile. Then in 
Sect. [6] we introduce the proximity effect strength parameter and 
develop additional models to disentangle various biases in the 
signal. In Sect. [7] we present the results for the mean line of 
sight proximity effect as determined from synthetic spectra. In 
Sect. |8] we discuss the proximity effect strength distribution of 
our spectra. Finally we summarise our findings in Sect. [9] 

2. Simulations 

2. 1 . Initial Realisation (IGM Model) 

We employ a DM simulation in a periodic boxQ] of 50 Mpc h~ l 
with 512 3 DM particles dvon Benda-Beckmann et ai] [2008). 
Using the PM-Tree code GADGET2 dSpringel et al]l2005l) our 
simulations yield a force resolution of 2 kpc hr l and a mass res- 
olution of m p = 7.75 x 10 7 Mq1i~ x . We fix the cosmological 
parameters to be in agreemen t with the 3rd year WMAP mea- 
surements (Sp ergel et al]|2007l) : The total mass density parame- 
ter is Q„,,o = 0.3, while the baryon mass density is Q^o = 0.04, 
and the vacuum energy is Qa = 0.7. The dimensionless Hubble 
constant is h - 0.7, and the power spectrum is normalised by the 
square root of the linear mass variance at 8 Mpc hr x , erg = 0.9. 

Simulating the Lya forest is a challenging task, since a 
large box is required to capture the largest modes that still in- 
fluence the mea n flux and its distribution dTvtler et al] 12009; 
iLidz et all l2010h . Furthermore high resolution is required to 
capture the properties of single absorbers correctly. Our simu- 
lations yield a reasonable repr esentation of the average statis- 
tics of the Lyman alpha forest dTheuns et al.ll 19981: IZhang et al. 



Distances are given as comoving distances unless otherwise stated. 
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Fig. 1. Cuts through the 25 Mpc /; 1 simulation box at redshift z — 4 for a QSO sitting in the void at the lower left corner. In the left 
panel the over-density field 1 + 6 is shown. The right panel provides the inverse of the over-ionisation fraction «Hi,Ly-<i/«Hi,prox due 
to a L Vll = 10 31 erg Hz ; s _1 QSO. Here «Hi,Ly -a and «Hi,prox are the neutral hydrogen fractions without and with the influence of 
the QSO. Clearly visible are the extended over-ionisation zone and shadowing effects. The scattered white dots arise from numerical 
noise. The axes are in comoving h Mpc. 



1 998> iMcDonald et~alll200lt IViel et alJl2002al) for low redshifts 



Ho wever recently it h as been shown by Bol ton & Becker] (2009) 
and lLidz et alj ( 120101) that to resolve the Lya forest properly at 
high redshifts, the mass resolution has to be increased by two 
orders of magnitudes for z = 5 and one order of magnitude for 
z = 4. Given the large box sizes needed, these resolution con- 
st raints cannot be fully m et with todays simulations. According 
to Bolton & Becker (2009), the error in our estimate of the mean 
flux at z = 4 is around 10%. For z = 3, the simulations will be 
marginally converged with mass resolution. 

We record the state of the simulation at three different red- 
shifts z = 4.9, 4, and 3. Using cloud in cell assignment we con- 
vert the DM particle distribution into a density and velocity field 
on a 400 3 regularly spaced grid. 

We then select three different environments from the highest 
redshift snapshot: The most massive halo with a mass of Mhaio = 
(2.6 x 10 12 ,4.0x 10 12 ,7.9x 10 12 , 1.0 x lO^M©^ 1 at redshifts 
z = (4.9,4,3,0), a random filament, and a random void. The 
last two environments are selected by visual inspection of the 
particle distribution and their position is tracked down to z = 3. 
These locations will be considered to host a QSO. 

To better characterise these environments, we estimate the 
volume-averaged overdensities 5s in r — 5 Mpc h spheres and 
find5 5 = (1.7, 1.8,2. 1) at redshifts z = (4.9, 4.0, 3.0) for the halo, 
(5 5 = (1.0,0.8, 1.1) for the filament, and S 5 = (0.8,0.7,0.7) for 
the void. Thus the halo resembles a locally overdense region, the 
void an locally underdense region, and the filament an average 
environment. Figure Q] illustrates a snapshot of the DM density 
field at redshift z — 4. 

In order to describe the baryonic component of the IGM, 
we assume a universe containing hydrogen only, whose density 
and velocity fields a re proportional to those of the dark matter 
dPetitiean et al.l 1 1 995h. In ord er to account for pressure effects 
on barvons. iHui et al.l ([1997) proposed to convolve the density 
field with a window functio n cutting off power below the Jeans 
length. In lViel et all d2002bl) differences of the gas and DM den- 



sities between hydrodynamical and DM only simulations are 
studied. For 6 ~ 3, baryons follow the DM distribution quite 
well, i.e. over most of the densities relevant for the Lyman-ff 
forest. However at higher DM densities, the corresponding gas 
density is lower due to the smoothing induced by gas pressure. 
We implicitly smooth our density field with the cloud in cell 
density assignment scheme and cell sizes of 125 h kpc. This 
is comparable to the Jeans length of w 150 hr x kpc at z — 3 and 
mean density, which scales as <S~'/ 2 (1 + z)~ 3 / 2 . 

The hydrogen density at position x is then given by 



« H (x) 



3 Hq 

8 7T G m„ 



(l+z) 3 (l+5(x)) 



(1) 



where G is the gravitational constant, Hq the Hubble constant, 
and trip the proton mass. The DM overdensity is 1 + 6(x) = 
p(x)/p, with DM density p(x). Finally the hydrogen velocity 
field v(x) is assumed to be equal to the DM one. 

2.2. Model and calibration of the intergalactic medium 

The thermal evolution of the IGM is mainly determined by the 
equilibrium of photoionisation heating and adiabatic cooling, re- 
sulting in a tight relation between the density and temperature 
of the cosmic gas. This relation i s known as the effective equa- 
tion of state dHui & Gnedinlfl997l) . and is typically expressed by 
T = 7o(l + 6) y ~ l . We apply the effective equation of state as an 
estimate of the baryonic density for overdensities 0.1 < S < 10 
relevant for the Lyman alpha forest. Larger densities concern 
collapsed regions which we approximate by assuming a cut-off 
temperature r c iit - n ff = T(6 = 10). Following the formalism of 
IHui et al.1 (1 19971) . we can compute an H i Lya absorption spec- 
trum from the density and velocity fields, once To, y, and the UV 
background photoionisation rate (Fuvb) are fixed. Thus through 
a match between simulated and observed H i absorption spectra 
properties, we can constrain these free parameters. Such a cali- 
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Fig. 3. Upper panel: Comparison of our choices for y (black 
points ) with observationally derived results by ISchave et alJ 
(2000) (gray points). Lower panel: C omparison of o ur mo del TO 
(black points) with observations by ISchave et al.l ((2000) (gray 
points). 



bration will be crucial in the following to provide realistic repre- 
sentations of the IGM for the radiative transfer calculations. 

To calibrate our spectra, we employ three observational 
constraints sorted b y increasing importance: (i) The obs erved 
equation of state dRicotti et al.l l2000t ISchave etaTI l2000h . (ii) 
the evolution of the UV background p hotoionisation rate 
(lHaardt&Madaull2001r. iBianchi et all l200lh . and (iii) the ob- 
served evolution of the effective op tical depth in the Lya forest 
( ISchave et alJl2003l;lKim e"tai]|2007l) . 

Our catalogues of synthetic spectra consist of 500 lines of 
sight at each redshift randomly drawn through the cosmological 
box. The spectra are binned to the typical resolution of UVES 
spectra of 6.7 km s and are further convolved with the UVES 
instrument profile. Random noise with a signal to noise (S/N) 
ratio 70 is added, the m ean S/N in the QSO sample used by 
iDaA'A grio et al. (2008a). An example of such a mock spectra 
is presented in Fig. [2] 

The model parameters of the effective equation of state Tg 
and y wer e chosen according to observations of iRicotti et alJ 
(2000) and lSchave etaTI (2000). To obtain values for z > 4, the 
observed To and y parameters have been extrapolated to higher 
redshifts (see Fig. [3]). We further constrained our model to yield 
an observed average effective optical depth r e ff(z) = - In (F(z)) 
where F is the transmitted flux, and the averaging is performed 
over the whole li ne of sight. For this we used recent observations 
from iKim et all (120071) . With these constraints we determined 
the UVB photoionisation rate Tuvb for our models. 

Our model parameters are presented in Table [TJ and are plot- 
ted in Fig.|4]in comparison with different literature results. Both 
the inferred evolution of the UV background and the effective 
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Fig. 4. Upper panel: The evolution of the UV background pho- 
to ionisation rat e in o ur three snapshots (black points) compared 
to iBolton etaTI (120051) (grey poin ts shifted by z = 0.05 fo r bet- 
ter visibility) and predictions by Haardt & Madau (2001) (grey 
line). Lower panel: The effective optical depth o f our models 
(black points) in comparison to measurements by ISchave et alJ 
d2003l) (g r ey po ints). The continuous line shows the fit by 
iKimet alJ J2007I) . 



Table 1. Model parameters of the semi-analytical model. 



z 


logTo 


7 


Tuvb 




log r cfI 


4.9 


4.23 


1.20 


0.6 x 10" 


12 


0.186 


4.0 


4.18 


1.15 


1.0 x 10- 


12 


-0.114 


3.0 


4.36 


1.10 


1.2 x 10 


12 


-0.474 



optica l de pth closely follow r ecent results by lHaardt & Madaul 
( 1200 lh andlBolton et all d2005l). and high resolution observations 
bv lSchave et al.l(l2003l) and lKim et alf(l2007l) . respectively. 



3. Radiative Transfer in Cosmological Simulations 

3.1. Method 

We now focus on the method implemented to solve the radia- 
tive transfer equation. In developing our Monte-Carlo radia- 
tive transfer code ( A-CRASH), we clo sely follow the approach 
first introduc e d by Ciardi et all (12001 1) and further extended by 
Mase lli et alJ (120031) . While their latest developmen t includes 
multi frequency photon packages (Ma selli et al.ll2.Q09b . we only 
implement the earlier formulation which uses monochromatic 
packages^. The advantage of the Monte-Carlo scheme is that it 
accounts for the emission of diffuse recombination photons, typ- 
ically neglected in codes relying on the "on-the-spot" approxi- 
mation (see Section l5~3"1 for further details). 

The idea behind a Monte-Carlo radiative transfer scheme is 
to bundle the radiation flux field into a discrete amount of radia- 
tion energy in form of photon packages having 



AE t = r 

•J',-1 



L\jy(t)dt 



(2) 



2 Our A-CRASH code is OpenMP parallel 
and publicly available under the GPL license at 
http : //sourcef orge . net/pro j ects/acrash/ 
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as energy content, sent out by a source with UV luminosity Ljjv 
between two time steps f,_i and % At each time step f, we allow 
the corresponding content of energy A£, to be evenly distributed 
among the number of photon packages emitted by the source. 
In our implementation of the Monte-Carlo scheme this number 
can be larger than one. Additionally each photon is characterised 
by a given frequency v following a given spectral energy distri- 
bution (SED). Thus, knowing the source SED and its UV lu- 
minosity we are able to infer the number of photons N 7t , per 
package. To guarantee a proper angular and spatial sampling of 
the radiative transfer equation, the photon packages are emitted 
in random directions. 

Once all the photon packages are produced by the sources, 
they are propagated through the computational domain. Every 
time a package crosses a cell, a certain amount of photons is 
absorbed. The absorption probability in the Z-th cell is 



P(t,) = 1 - e- T > 

with 77 being the optical depth in the Z-th cell 
Ti = cr Hl (v) n Hl j fiAx 



(3) 



(4) 



where cr H , is the hydrogen photo-ionisation cross-section, nm, / 
the neutral hydrogen number density, and fiAx the crossing path 
length of a ray through a cell with size Ax. We calculate the 
exact crossing length using the fast voxel traversal algorithm by 
lAmanatides & Wool dl987h . 

The number of absorbed photons in the cell is Na. / = 
N y . 1 P(ji) where N 7t 1 denotes the remaining number of photons 
in the package N 7t / = N y , /_i - Na, 1-1 arriving at the Z-th cell. In 
our simple case of pure hydrogen gas, we can write 



dXHn 

n H j = THtU )»Hi«e 

at 



■ «Hn(T)nHu«e +r H i«Hi 



(5) 



where ih„ = hhh/"h is the ionisation fraction, «e is the total hy- 
drogen density, ymiT) is the collisional ionisation rate, an u (T) 
the recombination rate, and Thi is the photoionisation rate de- 
rived from the number of absorbed photo ns Na, i- We adopt 
y H i(T) and a H n<7) as in lMaselli et ail ((20031) . We solve this stiff 
differential equation using a 4th order Runge-Kutta scheme. 

The extent of a time step is defined as the total simulation 
time t s divided by the total number of photon packages N p emit- 
ted by each source. The numerical resolution of the simulation is 
determined by calculating the mean number of packages cross- 
ing a cell 



N a 



N S N P 
N} 



(6) 



where N & is the number of sources and N c the number of cells 
in one box dimension, and f m i n is the smallest characteristic time 
scale of all the processes involved. In order to efficiently paral- 
lelise the scheme, the original rule of producing one photon per 
source per time step is dropped; a source is allowed to produce 
more than one package per time step as long as At <K f m j n . 

We have confirme d that our implemen tatio n passes the s im- 
ple tests described bv lMaselli etalj d2003l) (see lPartll d2007l) for 
further details). 

3.2. UV Background Field 

The UV background photoionisation rates given in Table Q] 
need to be modelled in the framework of the radiative trans- 
fer scheme. We assume the spectral shape of the UVB to be a 



power law v ab with at, = -1.3 dHuietalJ 1997b which is a bit 
lower than recent measurements bv lFechnereTafl (|2006) yield- 
ing ab = -1.99 + 0.34. From the spectral energy distribution of 
the UV background field and its intensity, photon packages need 
to be constructed. To obtain the number of photons in a back- 
ground package, the total energy content of the background field 
in the box is mapped to single phot on packages; we follow the 
method by iMaselli & Ferraral d2005l) . but for the hydrogen only 
case. 

To derive the energy content of a single background field 
photon package N y , the total amount of energy carried by the 
background field in the whole simulation box needs to be con- 
sidered. The change in the mean ionised hydrogen density «Hn 
by the UVB is 



An 



H11 



UVB 



«HiAf, 



(7) 



where Tuvb is the photo-ionisation rate of the background field. 
This is compared to the total amount of absorbed photons in the 
box N 7 (1 - exp(-Ar)), where At = cr Hl A«hi <^b is the total 
optical depth in the box of length d%. This corresponds to a mean 
change of the ionised hydrogen number density of 



An 



Hn 



N r (l -e- AT ) 
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(8) 



By equating Eqs.[7]and[8]and assuming an optically thin medium 
At <k 1 with At * (l - <?~ At ) we obtain the photon number 
content of a background photon package 



Ny = 



r UVB Af(Ax-AQ 2 

0"Hi 



(9) 



Note that this is only true, if the photon package is propagated 
over the distance of exactly one box length. 

Background photons are emitted isotropically from random 
cells in the box. Dense regions are allowed to shield them- 
selves from the background flux. Thus we emit background pho- 
tons only from cells bel ow a certain density threshold £uv- In 
Mas elli & Ferraral {2005) a threshold of <Suv = 60 was used, cor- 
responding to the density at the virial radius of collapsed haloes. 
We chose a lower threshold 6uv — 1 in order to ensure that 
mildly overdense regions have the possibility to shield them- 
selves from the UV background. 

3.3. QSO radiation 

As with the UV background, discrete point sources are charac- 
terised by a SED and luminosity. We only include one source 
of radiation other than the UVB. The point source represent- 
ing the QS O was chosen to f ollow the composite QSO spectra 
obtained by Tram mell et al.l d2007l) . The mean SED was con- 
structed from over 3000 spectra available in both the Galaxy 
Evolution Explorer (GALEX) Data Release 1 and the Sloan 
Digital Sky Survey (SDSS) Data Release 3. It covers a wide 
wavelength range of about 9000 > A > 300A. 

If a power law v a i is assumed for the U V part of the QSO 
spectrum, the data of Tramm ell et al.l ((2007) leads to a q « -2.5 
for A < 912A. Since the scatter in the data of Trammel l et al.l 
d2007h is large in the wavelength interval we are interested 
in, and other au thors find harder spectra dTelfer et alJ [2002; 
IScott et al] 120041) . we consider in addition a power law SED 
with a shallower slope of a q = -1.5. The upper energy limit at 
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Fig. 5. Lya absorption along a sight line near a z — 3 L Vll = 
10 32 ergHz"' s- 1 

QSO simulated with photon numbers N p q S0 = 1 x 10 s (solid line), 
1 x 10 9 (dashed line), 2x 10 9 (dotted line), and 4x 10 9 (dash dotted line). 
The wavelength is in Angstrom and no noise is added to the spectra 
for better visibility. The lower panel shows the relative deviation to the 
highest resolution run. 



A = 300A results in a underestimation of the hydrogen photoion- 
isation rate of 0.6% for a = -1.5 and only 0.2% for a = -2.5 
with respect to spectra extending to higher energies. 

The two QSO Lyman limit luminosities L Vll studied 
chosen to bracket the luminosity range of observed 

i et alj 
QSOs 



are cnosen to bracket tne luminosity rang e oi oDse 
QSOs dScott et al J 120001: iRollinde et ail 120051: rGuimaraes i 
120071: lDairAghoetalJl2008al) . We therefore chose the Q 



to have Lyman limit luminosities of 10 31 erg Hz 1 s 1 and 
10 32 erg Hz s _1 . In order to obtain the total UV luminosity for 
use with Eq. [2] the high energy part of the SED below 912A 
is scaled to the given Lyman limit luminosity and integrated 
from 912 > A > 300A. For the lTrammell et all d2007l) spec- 
tral template and L VLL = 10 31 erg Hz 1 s _1 , we obtain a total 
UV luminosity of L uv = 1.15 x 10 46 erg s" 
10 32 erg Hz" 1 s" 1 a luminosity of L uv = 1.15 x 10 47 = 



and for L VhL 
7 erg s _1 . 



3.4. Simulation Setup and Convergence 

Due to the high precision required to resolve the slight changes 
in the ionisation fractions governing the proximity effect, a large 
amount of photon packages had to be followed. To achieve this, 
sub-boxes of 25 Mpc hr x with 200 3 cells have been cut out of 
the whole box. To ensure the resulting proximity effect region to 
be traced to as large a radius as possible, the QSO source was 
located at the box origin at the cost of analysing only 1/8 of the 
full sphere around the QSO. 

A convergence study is carried out for one of our models 
with resulting mock spectra shown in Fig. [5] We test QSO sam 
plings JVp,Q S0 = 10 8 , 10 9 , 2 x 10 9 , 4 x 10 9 for a L VhL 
10 32 erg Hz" 1 s _1 QSO residing in a filament at z = 3. For this 
test the a = -2.5 spectrum has been used. Simulating the bright 
QSO at z — 3 provides the most challenging convergence test. 
The IGM is already highly ionised, and the influence of the 
quasar result in very small neutral hydrogen fractions and thus 
need high numerical resolution. The lower luminosity QSO will 
not alter the IGM as strongly and therefore resolving the ioni- 
sation fractions numerically is easier. The same is true at higher 
redshift, where the IGM is not yet as strongly ionised as at the 



low redshifts, and the resulting neutral fraction are larger due 
to the larger optical depth. The a = -2.5 spectrum is used for 
this test, since it is more challenging to properly sample the high 
energy tail compared to the harder a = - 1 .5 model. 

It is obvious, that 10 8 photons for the QSO are insuffi- 
cient, even if the mean number of photons crossing each cell is 
around 10. This number has been considered suffic ient to resolve 
Strom gren spheres in a homogeneous medium (Masel het alj 
120031) . The solutions for N p , QS o = 2 x 10 9 and 4 X 10 9 are sim- 
ilar, with relative deviations to the highest resolution run of at 
most 5 per cent. In order to achieve good angular sampling of 
the QSO environment and sufficiently large cell crossing num- 
bers we consider the solution for A^.qso = 2 x 10 9 as converged. 
To correctly sample the whole sphere, 1.6x 10 10 packages would 
be required, which is at the moment beyond the capabilities of 
our code. 

The Monte-Carlo radiative transfer method is able to follow 
photo-ionisation heating, however since this effect is already in- 
cluded in the equation of state, we keep the temperatures de- 
rived from the equation of state fixed throughout the simula- 
tion. This approach neglects any temperature fluctuations due 
to He ii photo-heating during He n reioni sation around z ~ 3 
(McOuin nlst al.ll200 9: Meiks in et aT]l2010h . The simulations are 
run up to t s = 2.5 x 10 8 yr and the ionisation fractions were av- 
eraged over different outputs at different time steps to reduce the 
Monte-Carlo noise. 



4. Line-of-Sight Proximity Effect 

In addition to the radiative transfer method for simulating the 
proximity effect, we also exploit a standard approach to account 
fo r the QSO radiation field. This me thod builds on th e work 
of iBaitlik et ail (Il988l also BDO) and iLiske & Willigerl (l200lh 
and we briefly summarise it here. Under the assumption that the 
IGM is in photoionisation equilibrium with the global UV radia- 
tion field, the amount of UV photons increases in the vicinity of 
bright QSOs and dominates over the UVB. This effect leads to 
a reduction in neutral hydrogen, and with it to an opacity deficit 
of the IGM within a few Mpc from the QSO. 

The influence of the QSO onto the surrounding medium 
leads to an effective optical depth in the Lycc forest 



Teff,QS0(z) = Teff.Lyo-fe) (1 + 0)(z)) 



1-/5 



(10) 



dLiske & Willigerl 1200 lh where T e ff,Lyo-(z) represents the evolu- 
tion of the effective optical depth in the Lya forest with redshift, 
and r e ff i Qso(z) is the effective optical depth including the alter- 
ations by the QSO radiation. Further /3 is the slope in the column 
density distribution and a>(z) is the ratio between the QSO and 
the UVB photoionisation rates. 

Then following the assumption of pure geometr i cal d ilution 
of the QSO radiation as proposed by Baitlik et alj d!988l) . we 
obtain 



= ^BDO(Z) = 



167T 2 J VLL (1 +Z)d L (Za,Z) 2 



(ID 



with z as the redshift of absorbers along the LOS such that z < Zq, 
and di{z q ,z) is the luminosity distance from the absorber to the 
QSO. L Vll is the Lyman limit luminosity of the QSO and / VLL 
the UVB flux at the Lyman limit. In observations, the emission 
redshift of the QSO z q is subject to uncertainties which propa- 
gates through to the u scale and increa se the uncertainti es in the 
determined UVB photoionisation rate dScott et al. 2000). In this 
work we consider the emission redshift to be perfectly known. 
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Fig. 6. Median over-ionisation profile as a function of the distance from the QSO (black line, grey area gives upper and lower 
qualities) hosted by a halo at z — 3.0 (left) and z = 4.9 (right side) for a QSO luminosity L YLL = 1 x 10 31 erg Hz -1 s _I . The analytical 
over-ionisation profile including the SED effects is given by the solid line without the grey area. The region affected by oversampling 
(r < 1 Mpc ) is marked with a dotted line. 



We note that this formula is derived assuming identical spec- 
tral indices of the QSO and the UVB. In our simulation we con- 
sider two different SEDs for the UVB and the QSO. Therefore 
in using Eq. QT| to estimate the proximity effect, we would in- 
troduce a bias. In u>(z) the two different spectral shapes can be 
accounted for by the ratio of the photo-ionisation rates of the 
QSO as a function of radius Fqso(>") and the photo-ionisation 
rate of the UVB, r UV B 



, . r QS o(r) 3 - a h 
«(z) = — = wbdo(z) ■ 



(12) 



UVB 



where at, and a q describe the UVB and the QSO SED, re- 
spectively, by power laws f v oc v° . Further it is assumed that 
<t Hi oc v -3 . Finally, the to scale is uniquely defined once the QSO 
Lyman limit flux, redshift, and the UVB flux are known. In the 
following analysis we will always use Eq. [T2] unless otherwise 
stated. 

This correction for differing SEDs has recently been 
used in measurement s of the UVB photoionisation rate 
dDairAglio et alJ |2009). However depending on the spectral 
shape of the QSO, omitting this correction can result in over- or 
underestimation of the background flux. As stated in Section [3~3l 
we use a QSO spectrum with a rather steep slope of a q = -2.5. 
Just using the original BDO formulation would result in an over- 
estimation of the UV background flux by almost 30% assuming 
ffb = -1.3. We will later check if this correction is able to model 
the SED effect by comparing our simulations with one using a 
shallower QSO spectral slope of a q = -1.5. 

As a side note, this spectral effect will alter the QSO's 
sphere of influence. We take as size of the proximity ef- 
fect the radius r„ = i where Tqso = Tuvb- Then r w= \ oc 
>"bdo,«=i -y/(3 - «b)/(3 - a q ), and the radius is in the case of our 
steep QSO spectra 13% smaller for equal L VLL . 



5. Radiative Transfer Effects 

5.1. The Over-ionisation Profile 

The over-ionisation profile of the proximity effect around the 
QSO decays radially due to geometric dilution. Additionally, 
differences in the spectral energy distribution between the UVB 



and the QSO influences the size of the proximity effect zone. 
Both factors are included in the o> scale. To determine how 
much the QSO decreases its surrounding neutral hydrogen den- 
sity in our simulations we determine the over-ionisation profile 
/«Hi,Prox, where HHiXy-o- gives the neutral hydro- 
gen fraction unaffected by the QSO radiation. Further «Hi,Prox 
denotes the neutral hydrogen density in case of including the 
additional QSO radiation. Assuming ionisation equilibrium, the 
over-ionisation profile is directly proportional to the u> scale. 

In Fig. [6] we show the median over-ionisation profile deter- 
mined on 200 lines of sight. These have been extracted from the 
simulation of a L Vll = 10 31 erg Hz s _1 QSO residing in the 
most massive halo at z = 3 and z = 4.9. For both redshifts the 
simulation data closely follows the analytical model described 
in Section|4]at radii r ^ 5 Mpc hr l at z — 3 and r ^ 2 Mpc hr l at 
z = 4.9. However, the faint QSO, the simulated over-ionisation 
profile starts to deviate at smaller radii, where the near host en- 
vironment leaves an imprint in the profile. In contrast, for the 
case of the bright QSO shown in the upper left panel of Fig. [8] 
the profile follows better the analytical model. We need to em- 
phasise that the spatial resolution of the simulations are limited 
to r — 0.125 Mpc hT 1 . Therefore the immediate vicinity of the 
QSO is badly resolved. Furthermore, cells up to 1 Mpc are sub- 
ject to oversampling in all our cases. Only beyond 1 Mpc cells 
are at most sampled once every time step. Therefore the solution 
in this part of the over-ionisation profile, which is marked in our 
plots, is unrealiable. 

5.2. Shadowing by Lyman Limit Systems 

We now focus on the over-ionisation field and discuss how op- 
tically thick regions in the intergalactic medium affect the prox- 
imity effect. A slice of the over-ionisation field is shown in the 
right panel of Fig. [Qfor the L Vll = 10 31 erg Hz 1 s" 1 QSO resid- 
ing in a void at z = 4. Again the smooth over-ionisation profile is 
seen as expected. However in this smooth transition zone, large 
shadow cones originating at dense regions are present. The op- 
tical depth of these clouds is large enough to absorb most of 
the ionising photons produced by the QSO. Hence the ionisa- 
tion state of the medium behind such an absorber stays at the 
unaltered UVB level. 
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Fig. 7. Mock spectra synthesized from the z — 4 simulation with 
the QSO in a void. The solid line gives the Lya forest spectrum 
without the influence of the QSO, the dashed line the spectrum 
from the radiative transfer simulation of the proximity effect for 
a QSO Lyman limit luminosity of L Vll = 10 31 erg Hz s _1 . The 
grey line gives results from a semi-analytical model of the prox- 
imity effect showing strong deviations from the radiative transfer 
results behind the strong absorber at the middle of the spectrum. 
No noise is added to the spectra for better visibility. 



We determined the hydrogen column density of the shadow- 
ing regions in Fig.Q] They lie in the range of 1 x 10 17 < Nm < 
5 x 10 18 cirr 2 . The regions causing the shadows therefore rep- 
resent Lyman limit systems. Whether such systems are able to 
absorb a large amount of QSO photons depends on the one hand 
on their density and on the other on the QSO flux reaching them. 
Absorbers further away from the source will receive less flux 
from the QSO due to geometric dilution. Therefore they are more 
likely to stay optically thick. For systems nearer to the QSO, the 
flux field is larger and a higher density is needed for them to 
remain optically thick. Thus the number of shadows increases 
with increasing distance to the source. Moreover a highly lumi- 
nous QSO will produce a larger amount of ionising photons and 
a higher number of absorbers are rendered optically thin. Hence 
less shadows are present in the over-ionisation field of a highly 
luminous QSO. 

By extracting a spectrum through such an absorber (see Fig. 
UJ) the influence of the Lyman limit system on the proximity ef- 
fect signature can be seen. Behind the strong absorber on the left 
hand side, the spectrum follows the unaffected Lya forest spec- 
trum. Measurements of the proximity effect behind such a sys- 
tem do not show any QSO influence. This shadowing is not de- 
scribed with a semi-analytical model where the optical depth in 
the forest is altered proportional to 1 + (jl>(z). We will discuss and 
make use of this semi-analytical model in Section|6] By looking 
at the semi-analytic model spectrum shown in Fig. [7] it is clear 
that in this model the proximity effect extends beyond the Lyman 
limit system. 

5.3. Diffuse Recombination Radiation 

We will now discuss the role played by photons produced by 
recombining electrons with the help of the over-ionisation pro- 
file, since our radiative transfer simulations give the opportunity 
to study this diffuse component. With the Monte-Carlo scheme 
we can directly model the radiation field produced by recom- 
bination, which in most other schemes is not treated self con- 
sistently. Diffuse radiation produced by recombining electrons 
plays an important role in radiative transfer problems if the 
medium is optically thin. Whenever photons originating in re- 



combination events are able to travel over large distances (i.e. 
larger than one cell size), energy is redistributed on the scale of 
the photon's mean free path and the ionisation state is altered 
over these distances. The importance of diffuse radiation in the 
simplest of radiati ve transfer prob lems, the Stromgren sphere, 
has been shown in iRitzerveldl (120051) where the outer 30% of the 
Stromgren radiu s are found to be dominat ed by recombination 
photons. Further lAubert & Tevssi er (2008) used their radiative 
transfer code to discuss how much the on-the-spot approxima- 
tion, a popular assumption that treats recombination photons as 
a local process, alters results in widely used test cases. 

Recombination photons to the hydrogen ground level pos- 
sess enough energy to ionise other hydrogen atoms. If such a 
photon is absorbed by a nearby neutral atom, the on-the-spot 
approximation can be applied where it is assumed that all re- 
combination photons are absorbed by neighbouring atoms. This 
is true if the medium is dense and the optical depth high enough 
that a local equilibrium between recombination and absorption 
is reached. If however the optical depth is small, the mean free 
path is larger than the simulation cell size and the recombination 
photons are absorbed far away from their origin. In this case, the 
on-the-spot approximation breaks down. 

We expect recombination radiation to play an important role 
in the proximity effect, since the optical depths in the Lya forest 
are low. To identify the region where recombination processes 
influence the solution, we run one simulation without emitting 
recombination photons. However we still use case A recombi- 
nation rates and call this the 'no diffusion case'. Furthermore 
we implemented the on-the-spot approximation using case B 
recombination rates. As before we do not follow any recom- 
bination photons. We expect to see differences at lower red- 
shift, where the mean free paths for ionising photons are large. 
Therefore we consider the snapshot at redshift z — 3 with the 
L Vll = 10 32 erg Hz 's 1 QSO residing in a halo for this experi- 
ment. 

The resulting median over-ionisation profiles are presented 
in Fig. [8] We also show the over-ionisation distribution at two 
radii to better compare the different runs. At radii r > 1 Mpc h~ l 
we obtain profiles that well reproduce the theoretical estimates. 
We see some noisy fluctuations in the profiles. Due to the larger 
ionisation fractions for the bright QSO, the ionisation fractions 
themselves are more susceptible to Monte-Carlo noise. Since the 
neutral hydrogen fractions in the direct vicinity of the QSO are 
2 dex lower for the bright QSO than for the faint one, higher 
precision is needed to better evaluate the extremely low neu- 
tral fractions. However further increasing the number of photon 
packages is at the moment beyond the capabilities of our code. 

Comparing the no diffusion run with the full radiative trans- 
fer solution shows that recombination greatly contributes to the 
outer parts of the proximity effect region. The no diffusion solu- 
tion does not gradually go over to unity, but keeps on decaying 
fast until the QSO's influence vanishes. At r — 20 Mpc h~ l the 
difference between the result with and without the diffuse recom- 
bination field is 23%. This means, that recombination radiation 
contributes to the photoionisation rate by 23% at this radius, if 
one assumes that the ionisation fraction is directly proportional 
to the photoionisation rate. 

The on-the-spot solution on the other hand over predicts the 
extent of the proximity effect zone by about a factor of 2. This 
means that keeping the energy produced by recombining elec- 
trons locally confined over predicts the amount of ionised hy- 
drogen. As a result, the QSO's sphere of influence is overesti- 
mated. Further the on-the-spot solution shows larger dispersion 
between the lines of sight than the radiative transfer simulation 
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Fig. 8. Median over-ionisation profile as a function of the distance from the QSO hosted by a halo as in Fig. [6] We show the influence 
of the diffuse field on the H n region for L Vll = 10 32 erg Hz ; s _1 at z = 3.0. The upper left plot shows the full radiative transfer 
simulation, the right plot the one without diffusive radiation, and the lower left plot results using the on-the-spot approximation. 
The analytical over-ionisation profile including the SED effects is given by the solid line. The region affected by oversampling 
(r < 1 Mpc ) is marked with a dotted line. The lower right histograms gives the distribution of the over-ionisation fraction at radii 
r — 10 Mpc h' 1 (left) and r — 20 Mpc hr x (right). The solid lines give radiation transfer results, the dashed lines shows results 
omitting diffusion, and the dot-dashed distribution illustrates results obtained with the on-the-spot approximation. The solid tick 
gives the analytical solution. 



(see Fig. [8] lower right panel). We therefore conclude that the 
on-the-spot approximation is not sufficient when modelling the 
proximity effect. 

6. Radiative Transfer on the Line-of-Sight Proximity 
Effect 

6. 1 . Strength Parameter 

In this Section we discuss the imprint of radiative transfer on the 
line-of-sight proximity effect in Lya spectra as analysed in ob- 
servations. To this aim we study a sample of 500 lines of sight 
originating at the QSO position and having randomly drawn di- 
rections. We then measure the proximity effect signature. 

First an appropriate w-scale is constructed for each QSO ap- 
plying Eq.[T2lincluding the SED correction term. In Section lTTSl 
we discuss how well the correction term accounts for the SED ef- 
fect. Given the w-scale we then determine the mean transmission 
in bins of A log to = 0.5 for each line of sight. The mean proxim- 
ity effect signal is calculated by averaging the mean transmission 
in each log to bin over the QSO lines of sight. We determine the 
mean optical depth in each bin. The relative change introduced 



in the optical depth by the QSO is evaluated by computing the 
normalised optical depth £ defined as 

£ = l55^, (13) 

' eff , Lya 

where T e ff,Qso is the effective optical depth in a logw bin and 
T e ff,Lyff the effective optical depth in the Lya forest without the 
influence of the QSO. For Alogw < 1 we find the resulting 
normalised optical depths to be weakly dependent on the chosen 
bin size. 

To quantify any difference in the proximity effect signature 
to the expecte d one given by Eq. [T2l we adopt the 'strength pa- 
rameter' as in lDairAglio et all d2008bllal) . The strength parame- 
ter a is defined as 

f = (l + -) (14) 

where to is defined by Eq. [T2] As we know the UVB photoion- 
isation rates in our models, we use the values given in Table Q] 
to define the reference to. Then values a > 1 or a < 1 describe 
weaker or stronger proximity effects, respectively. The SED cor- 
rection term in Eq.[T2]can also be interpreted in the framework 
of the proximity strength parameter. 
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It is also possible to determine the strength parameter for 
each line of sight individually. By obtaining the proximity effect 
strength for each line of sight, the strength parameter distribu- 
tion function can be constructed. Studying its properties gives 
further i nsight into the p r oximity effect signature and has been 
used bv lDaU'Aglioet all d2008bl) to derive the UVB photoioni- 
sation rate. How well this method performs in the context of the 
present study will be discussed in Section[8] 



6.2. Additional Models 

In addition to the lines of sight obtained from the radiative trans- 
fer simulation, we explored two other models allowing us to 
characterise the imprint of the cosmological density distribu- 
tion and radiative transfer in the analysis of the proximity effect. 
These complementary analyses will be used to obtain spectra 
that can be directly compared to our radiative transfer results. 



Semi-Analytical Model (SAM) of the Proximity Effect: From 
the radiative transfer simulation of only the UV background, we 
know the ionisation state of hydrogen in the box without the 
influence of the QSO. Then we can semi-analytically introduce 
the proximity effect along selected lines of sight from the QSO 
by decreasing the neutral fractions by the factor (1 + to). We 
choose the same 500 lines of sight as in the full radiative transfer 
analysis. 

This model includes the imprint of the fluctuating density 
field on the spectra. Since the semi-analytical model only in- 
cludes geometric dilution and the influence of the SED, any dif- 
ferences to the results of the full radiative transfer simulations 
reveals additional radiative transfer effects. 



Random Absorber Model (RAM): In our radiative transfer sim- 
ulations, we can only cover the proximity effect signature up to 
log to > 0.5 for the faint QSO and to log to > 1 .0 for the strong 
one, due to the finite size of our simulation box. In order to over- 
come this limitation, we employ a simple Monte Carlo method 
to generate Lye mock spectr a used in observations to study sys- 
temat ic effects in the data dFechner et al.l l2004; Wo rseck et al.l 
120071) . The Lyor forest is randomly populated with absorption 
lines following observationally derived statistical properties. The 
constraints used are: 1) The line number density distribution, 
approximated by a power law of the form dn/dz oc (1 + z) y 
where y = 2.65 (iKim et al.l 120071) . 2) the column density dis- 
tribution, given b y f(Nni) <* where the slope is /? 1.5 
dKim et al.ll2001l) . and 3) the Doppler parameter distribution, 
given by dn/db oc b~ 5 exp [-b 4 //? 4 ], where b a - 24 km/s 
(IKim etal.ll 2001). For a detailed discussion of the method used, 



see 



DairAglio etaj] (l2008allbh . 



The random absorber model does not take clustering of ab- 
sorption lines due to large scale fluctuations in the density distri- 
bution into account . Therefore we use this model to determine 
any possible bias in the analysis method and to infer the effect of 
absorber clustering. Furthermore, since we can generate spectra 
with large wavelength coverage as in observated spectra, we can 
quantify the effect introduced by the truncated coverage of the to 
scale in the spectra from the radiative transfer simulation. With 
this simple model it is possible to study any systematics present 
in the analysis, such as the influence of the A log to bin size on 
the resulting proximity effect. 



7. Mean Proximity Effect 

Now we discuss the results of the mean proximity effect as mea- 
sured in the synthesized spectra. 



7. 1 . Mean Normalised Optical Depth 

The main results are presented in Fig.[9]and Fig. [12] In both plots 
we show the analytical model according to Eq.[10]and Eq.[T3lfor 
comparison. The model predicts that at high to values, the QSO 
shows the strongest influence and the normalised optical depth is 
very low (at the QSO where to — oo, the normalised optical depth 
is £ = 0). Then between 2 > log to > -1 the influence of the 
QSO gradually declines due to geometric dilution and reaches 
£ = 1 for small to values. There the Lya forest is unaffected by 
the QSO. 

The normalised optical depth profile derived from the trans- 
fer simulations closely resembles the functional form of the an- 
alytical expectations. However a shift of the profile is present 
which we quantify using the proximity strength parameter. In 
observations, the strength parameter is used to measure the UVB 
photoionisation rate. We therefore faithfully apply the same 
method to the simulation results. However, since we know the 
UVB in our models, we can directly compare the photoionisation 
rate determined using the strength parameter with the model in- 
put. The strength parameter is determined by fitting Eq.[l4]to the 
data points using the Levenberg-Marquardt algorithm and a least 
absolute deviation estimator weighted with the error of the mean. 
In order to obtain a reasonable fit, data points deviating strongly 
from the profile above log to > 2 have been excluded. This also 
excludes the unreliable region within 1 Mpc of the QSO in the 
radiation transfer simulation, which translates to an log to > 2.4 
for the faint QSO and log to > 3.4 for the bright one. We provide 
the resulting strength parameters in Tableland include the fitted 
profiles in all plots of the mean proximity effect. The origin of 
the shift is due to the large scale environment around the QSO 
hosts which we discuss in more details in Section l731 

For each data point we show two error bars, the standard de- 
viation of the sample and the 2<x error of the mean. The sample 
standard deviation gives the variance between the 500 lines of 
sight. These variances are very large and extend to values above 
£ > 1 . Values of £ > 1 arise when saturated absorption systems 
dominate the optical depth in a log u> bin and raise the optical 
depth over the mean effective optical depth in the forest. A de- 
tailed illustration is given in the Appendix. Since the real space 
distance covered by the log to bins increases with decreasing to, 
there the contribution of saturated systems and the variance is 
reduced. 

The error of the mean is on the other hand small due to 
our sample size of 500 lines of sight, thus the mean profiles 
are well determined. However the mean profile does not strictly 
follow the smooth analytical model, but shows strong fluctua- 
tions. These fluctuations are enhanced by the fact that we can 
only analyse the QSO in one octant. Further we are focusing our 
analysis just on a single object, making the signal sensitive to 
the surrounding density distribution. If multiple halos would be 
combined to one mean profile, the fluctuations would diminish 
and the profile would follow the smooth analytical profile more 
closely. Since the data points deviate from the analytical profile, 
the determination of the strength parameter is quite uncertain. To 
estimate a formal error in the strength parameter we applied the 
Jackknife method. The error values are provided in Table [2] 

We checked if Monte-Carlo noise of the radiation transport 
simulation contributes to the signal's variance by reducing the 
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Fig. 9. The normalised optical depth for the L VLL = 
10 31 erg Hz 1 s _1 QSO with the standard SED and residing in 
a halo at z = 3. Continuous lines show the fitted profiles. The 
triangles and dash-dotted line give the radiation transfer proxim- 
ity effect without the use of the SED correction term. The cir- 
cles and dotted line give the radiation transfer results including 
the SED correction. The QSO with the hard a q = -1.5 SED is 
given by the squares and dashed line. The original BDO model 
is shown as the solid line. The dotted error bars show the sam- 
ple standard deviation. The solid errors show the 2cr error of the 
mean. 



QSO's photon packages production by half. This does not con- 
tribute noticeably to the variance in the signal. The large scatter 
in the proximity effect signal therefore has its origins in the dis- 
tribution of absorption systems along the line of sight. It is a 
direct imprint of the cosmic density inhomogeneities. The prob- 
ability distribution of the normalised optical depths in each log co 
bin is approximately log-normal as shown in Fig. IA.2l 

7.2. Performance of the SED Correction Term 

With our radiative transfer simulations it is possible to see, 
whether the SED correction term in Eq.[T2]describes the effect of 
a QSO and UV background which have different spectral shapes. 
The SED adds to the mentioned shift in the normalised optical 
depth, already introduced by the large scale density distribu- 
tion around the QSO hosts. In order to single out the SED influ- 
ence, we re-simulate the z — 3 faint QSO residing in a halo with 
a of q = —1.5 power-law SED. With the background's SED slope 
of ab = -1.28, we expect only a sl ight shift due to t he diff erent 
slopes of A log a = 0.02. The softer frrammell et al.1 d2007h SED 
used in the rest of the paper yields a shift of A log a = 0.11. If the 
SED correction term performs as intended, the corrected prox- 
imity effect profile of the soft SED QSO would approximately 
lie on top of the one with the hard SED. 

In Fig. [9] we present the results of this experiment, where 
we plot the proximity effect profile of the a q = -1.5 QSO, 
and SED corrected and uncorrected profile for the QSO with 
a q « -2.5. The data points of the corrected a = -2.5 profile lie 
almost on the a q - -1.5 profile as expected. A similar picture 
emerges by looking at the strength parameters. The uncorrected 
soft SED yields a strength parameter of log a = 0.84 + 0.06. By 
including the correction term we obtain logo = 0.70 + 0.12. 
For the a q = —1.5 SED we obtain a strength parameter of 
logo = 0.65 + 0.05. The different values of the SED corrected 
and the soft QSO results are consistent within the error bars. 
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Fig. 10. The mean overdensity profiles of the halo (solid line), 
filament (dashed line), and void (dash-dotted line) environment 
at redshift z — 3 as a function of comoving radius. The mean 
profiles are shown by thick lines. The thin lines give the Icr fluc- 
tuations in the density field. An excess in overdensity is seen 
for all environments between 3 < r c < 10 Mpc h~ l , where the 
profile stays above the average overdensity. 



7.3. influence of the Large Scale Environment 

The environmental density around the QSO has been shown 
to strongly influence the proximity effect. By studying ran- 
dom DM halos in a mass range of 1. 4 xIO'^q/i" 1 up to 
9 x lO l2 M Q h-\ Faucher-Giguere eTall (l2008h have found an 
enhancement in the mean overdensity profiles up to a comov- 
ing radius of (10 - 15) Mpc hT 1 . Only at larger radii the average 
overdensity profiles goes over to the cosmic mean. The overden- 
sity profiles fluctuate strongly between different halos. Since an 
enhanced density directly translates into an enhancement of the 
Lya forest optical depths, the proximity effect measurement is 
biased accordingly. 

We have determined the mean overdensity profiles together 
with the fluctuations around the mean for our three QSO host 
environments. The resulting mean profiles are shown in Fig.fTOl 
for z = 3. Our density profiles show two different influences. 

First, the local environment around the QSO host leaves an 
imprint on the profiles up to a radius of r ~ 2 Mpc h~ l . This 
is seen in the halo case where the massive host halo is re- 
sponsible for a strong overdensity which steadily declines up to 
r * 2 Mpc h~ x . The contrary is seen in the void case, where up 
to r » 3 Mpc h~ l the underdense region of the void is seen. At 
larger radii, the density increases quite strongly. 

All our environments show an overdensity between r » 
3 Mpc hT l up to a radius of around r ~ 15 Mpc hr l . In this re- 
gion the density lies above the cosmic mean. We refer to this 
phenomenon as a large scale overdensity. Further away from the 
QSO, the cosmic mean density is reached. The same behaviour 
is seen at higher redshifts. There, however, the amplitude of the 
large scale overdensity is lower. The presence of the large scale 
overdensity in the case of the most massive halo is a direct conse- 
quence of the fact that the most massive halo forms where there 
is a large overdensity. However around the filament and the void 
such a large scale overdensity is a random selection effect. 

An excess in density translates into an excess of Lyo- op- 
tical depth in this region, when compared to the mean optical 
depth. Thus the normalised optical depth £ increases and leads 
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Fig. 11. Comparison of the proximity effect at redshift z — 3 of a) the RAM (triangles and dotted line), b) the semi-analytical model 
(circles and dash dotted line), and c) the results from the radiative transfer simulation (squares and dashed line). The solid line shows 
the BDO reference analytical model. The dashed vertical line marks the start of the unreliability region of log co > 2.4 for the faint 
QSO and log co > 3.4 for the bright QSO. Results are shown for a QSO sitting in a filament. The left panel is for the low luminosity 
source, and the right one for the high luminous one. We only show the 2cr errors of the mean. 



to a biased proximity effect. The immediate environment around 
the quasar goes out to logw * 1.5 for the faint source and 
to logw « 2.5 for the bright one. At smaller logo; the large 
scale overdensity influences the proximity effect signal up to 
log co « 0.3 for the faint and to log co » 1 .3 for the bright QSO. 

Motivated by these considerations we study the influence of 
the cosmic density fluctuations on the proximity effect profile. 
To understand whether radiative transfer effects or the large scale 
overdensity are responsible for the shift in the normalised optical 
depth shown in Fig. [9] we make use of the reference models 
introduced in Section 16721 



We show the resulting proximity effect profiles for the refer- 
ence models of the faint and bright QSO residing in the filament 
at z = 3 in Fig.Qj] We use the RAM to infer the bias introduced 
by the analysis method itself due to the finite size of our spec- 
tra. Additionally the SAM allows us to determine the influence 
of the fluctuating density field by omitting any possible radia- 
tive transfer effect. On top of that we show the radiative transfer 
results. 

The data points of the RAM follow the standard BDO profile 
well for both luminosities. However a slight shift is seen in the 
fitted profiles. The strength parameters are log a - 0.06 + 0.05 
in the case of the weak QSO and logo = 0.05 + 0.10 for the 
strong QSO. Within the error bars, this is consistent with the 
analytical predictions. A small shift in the signal has been found 
to be intrinsic to th e analysis method due to the asymmetry in 
the ^-distribution dDalF Aglio & Gnedinll2009l) . 

As for the radiative transfer results, the SAM shows a strong 
shift towards the QSO. This is a clear signal, that the cosmo- 
logical density distribution is responsible for the large shift. 
Comparing the SAM with the radiative transfer results reveals 
no strong difference. The ratio between the strength parameter of 
the two models log(flRT/flsAivi) is for the faint QSO -0.1 and for 
the bright one 0.1. Considering the uncertainties in the strength 
parameters, the two model yield identical results. We provide 
the ratios of the strength parameters between the two models for 
all the studied environments and redshifts in Table [2] The ratios 
show a strong scatter between the different models. However due 
to the uncertainties in the strength parameter, we attribute these 
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Fig. 13. The combined proximity effect signature of the SAM 
applied to 500 randomly selected lines of sight with random ori- 
gin. A L VLL = 10 31 erg Hz s _1 QSO is assumed. The continu- 
ous line shows the fit of the analytical profile to the data. 



differences to the inability of the analytical model to account for 
the large scale density distribution. 

To strengthen this conclusion, we have performed one addi- 
tional test, where we applied the SAM to 500 randomly selected 
lines of sight with random origins and directions in the box and 
not just one specific origin. The resulting proximity effect pro- 
file is shown in Fig. [13] The profile now smoothly follows the 
analytical model and does not show any fluctuations around it. 
The strength parameter of log a = 0.1 1 + 0.10 is as well consis- 
tent with the SED corrected BDO model. The analytical model 
is thus a good description for a sample of lines of sight with 
random origins and directions. 

Thus the large shift in our simulations is solely due to the 
presence of large scale overdensities around our sources. Further 
we cannot identify significant radiative transfer effect in addi- 
tion to the one modelled with the SED correction term, since the 
SAM model gives similar results to the radiative transfer simu- 
lations. 
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Fig. 12. The combined proximity effect signature of all 500 lines of sight analysed with the BDO model showing the dependence 



on environment, luminosity and redshift. The left column shows the case of a faint L VLL = 10 erg Hz s QSO, the right one 
of a luminous L VLL = 10 32 erg Hz 1 s _1 QSO. From top to bottom redshifts are 4.9, 4, and 3. The solid lines show the geometrical 
dilution model, the other lines are best fits to the measurements with parameters given in Table [2] The circles and dotted lines 
show the halo cases, the squares and dashed lines the filament case and triangles and dash-dotted lines the void case. Symbols are 
horizontally shifted by +0.05 for better visibility. The dotted error bars show the sample lcr standard deviation. The solid errors 
give the 2<x error of the mean. The dashed vertical line marks the start of the unreliability region of log ai > 2.4 for the faint QSO 
and log io > 3.4 for the bright QSO. 



7.4. influence of the Host Environment 

We present the resulting proximity effect profiles for the halo, 
filament, and void QSO host environments as a function of red- 
shift in Fig.[T2] The corresponding strength parameters are given 



in Table [2] In the analysis the unreliability region of log a> > 2.4 
for the faint QSO and log oj > 3.4 for the bright QSO have been 
exculded. As we have seen in the previous Section, the local 
host environment of the QSO extends up to a comoving radius 
of r * 2 Mpc hT 1 . Hence, we expect the local environment to 
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Table 2. Parameters of the log a fits to the combined proximity 
profile for the semi-analytical model (SAM) and the radiative 
transfer simulation (RT). We give lcr errors. 
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affect the proximity effect profiles up to log w « 1 .5 for the faint 
QSO and to log a> w 2.5 for the bright one. 

The effect of the local environment is significant in the halo 
case. Except for the z = 4.9 bright QSO, the ^-profile rises 
strongly when approaching the direct vicinity of the QSO. The f- 
values start to rise to such an extent, that they lie out of the plot- 
ted region and deviate strongly from the analytical model. In ad- 
dition, the fluctuations in the normalised optical depth strongly 
increase. This departure from the analytic profile turns out to 
pose a problem in the determination of the strength parameter. 
In order to obtain a reasonable fit, the local environment needs 
to be excluded. 

For the filament and void environment however, such a de- 
parture is not seen as expected. Only in the profile of the z = 3 
faint QSO in the filament environment there is an enhancement 
at log to m 2.5. In the filament and void environments, the scatter 
in the normalised optical depth decreases. The lack of density 
enhancements in the immediate surroundings of the host thus 
reduces the scatter. The influence of the local environment also 
decreases with increasing luminosity, which is not only due to 
the reduced opacity, but as well due to the larger spatial cover- 
age of each log a> bin. 

Beyond the local environment's region, the ^-values return 
to the smooth profile expected from the analytic model. At 
z = 4.9 the different environments show identical profiles. With 
decreasing redshift however, the profiles in the different envi- 
ronments start to deviate from the analytical profile without any 
clear trend. These deviations are also seen in the corresponding 
strength parameters. For example the faint halo and void case at 
z = 4 coincide. Though for the bright QSO, the strength param- 
eters of the halo and the void do not agree anymore. Now the 
halo and the filament coincide within the error bars. Due to the 
limitations of the analytical model to account for the large scale 
overdensity present in the simulation data, the fitted strength pa- 
rameters show large unsystematic fluctuations. In Section [5] we 
have seen, that the local environment does not cause deviation 
from the geometric dilution model beyond 2 Mpc h~ . The ove- 
rionisation profile was not found to be influenced by the halo, 
filament, or void environment at large radii. Therefore we con- 
clude that the local environment does not influence the proxim- 
ity effect profile at radii larger than the local host environment. 
Furthermore the fluctuations in the strength parameter do not 
correlate with the source environment and solely arise from the 
fact, that the analytical model does not take the cosmological 
density distribution into account. 



7.5. Redshift Evolution 

Studying the strength parameters, no apparent evolution with 
redshift can be seen, if the fitting errors are taken into account. 
However the uncertainties of the strength parameter decrease 
with increasing redshift. Only in the faint halo case there is a 
hint of a decrease in proximity strength with increasing redshift. 

By looking at the ^-profiles in Fig.[T2]a clear decrease in the 
fluctuations around the mean profile can be seen with increasing 
redshift. This is due to the fact, that the density fluctuations of 
the dark matter are not as developed yet at high redshifts than 
at lower ones. A minor effect may result from the decreasing 
transmission of the Lya forest at higher redshift, i.e. the narrow 
range of gas densities probed by the Lyo- forest. The analysis of 
the large-scale overdensity as discussed in Section 17.31 demon- 
strates however the dominance of the evolution of the density 
field. There is a slight tendency of the profile to show a stronger 
shift at low redshifts than at z — 4.9. This is especially prominent 
in the case of the L VLL = 10 32 erg Hz 1 s -1 QSO. 

7.6. Luminosity Dependence 

The environmental bias of the proximity effect signature depends 
on the QSO luminosity. For all our redshifts, the shift in the pro- 
file is smaller for the bright QSO than for the weak one. Also the 
variance of the signal is reduced with increasing QSO luminos- 
ity. This is best seen for redshifts z = 4.9 and z = 3. 

The cause of the luminosity dependence lies in the fact that 
the u> scale is a function of the QSO luminosity. Therefore an 
increase in the QSO luminosity translates into larger distances 
covered by each log a> bin. This dampens the influence of the 
fluctuating density field and thus the variance in £ . Hence lumi- 
nous QSOs constrain the proximity effect signature better. 

8. Proximity Effect Strength Distribution 

Now we discuss our results on the distribution of the strength 
parameter determined on individual lines of sight. The asymme- 
try of the ^-distribution and the large scale overdensity biases 
the mean proximity effect profiles, as we have discussed in the 
previous Section. To avoid this bias, DalFAgli o et al.l {2008a) 
used the proximity strength distribution along individual lines of 
sight to measure the UVB flux, and proposed the peak of the dis- 
tribution to provide a measure of the UVB. We test the proposed 
method with our simulations. 

The proximity effect strength distribution is obtained by de- 
termining the ^-profile for each individual line of sight and by 
fitting Eq.[14]to the profile. In contrast to the last Section we use 
all available £ data points (except the ones inside the numerically 
unreliable region) in the fitting procedure. Whenever the simu- 
lated data deviates strongly from the analytical model due to the 
large variance of the ^-values, it is not possible to determine a 
strength parameter. This problem is significant in the case of the 
faint QSO residing in a halo, where the proximity effect profile 
deviates from the analytical one at log co > 1.5. There the fitting 
procedure failed in most cases and we cannot derive a strength 
parameter distribution. 

In Fig. [14]we present the proximity effect strength distribu- 
tion from our simulations. Only those lines of sight contribute 
to the distribution, for which a strength parameter could be de- 
termined, however we always normalise the distributions P(a) to 
unity. We employ the comparison models to investigate the var- 
ious effects that alter the strength distribution. In order to com- 
pare with the results obtained from the mean proximity effect 
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Fig. 14. The distribution of the strength parameter log a of the 
proximity effect in the different lines of sight for the radia- 
tive transfer simulations (solid histogram), compared with RAM 
(dotted line) and the SAM (solid line). The mean a-values from 
the radiative transfer simulations Tabled and their uncertainties 
are shown as vertical lines on top of each panel. The upper two 
panels show the low luminosity case. The lower three panels give 
results for the luminous QSO. From left to right the redshifts are 
4.9, 4, and 3. The vertical dashed line marks logo = 0.0. 



profile, we marked the strength parameters of the mean profiles 
as given in Table [2] 

All our models show a strength distribution that is clearly 
peaked. Furthermore the distribution reveals a large scatter of 
cr(loga) = 0.3, 0.4, 0.7 at redshifts z = 4.9, 4, 3. This scatter 
is much larger than the variance of the log a values of the mean 
analysis. All the models show a similar strength distribution as 
the respective ones of the RAM. In the case of the faint QSO, 
they all peak around log a = 0. For the bright QSO, the peak is 
somewhat shifted to negative log a = -0.1. This shift is as well 
present in the RAM. Therefore we suspect that it is due to the 
small a) scale coverage for the bright QSO and the corresponding 
uncertain fits. It follows that the peak of the strength distribution 
can be used to get unbiased estimates of the proximity effect. 



Surprisingly the strength distribution is thus not sensitive to 
the large scale overdensity environment discussed in the mean 
analysis. In the averaging process for deriving the mean profile, 
we are biased by the spherically averaged overdensities shown 
in Fig. [10] By fitting the strength parameter individually, less 
weight is given to single log u bin s deviating from t he exp ected 
profile. This confirms the claim in Dall'Agl io et al.l (l2008al) that 
using the proximity effect along individual lines of sight leads to 
a stable unbiased estimate of the proximity effect and thus of the 
UVB. 

As in the mean analysis, we can investigate the influence of 
the local environment, the quasar redshift, and the quasar lumi- 
nosity. The shape of the strength parameter distribution is not 
influenced by the QSO's local environment. Whether the QSO 
resides in a halo, a filament or a void, the maximum, median and 
the width of the distribution is not affected. Only for the faint 
quasar, the halo environment strongly influences the proximity 
effect profile so that in many cases no reasonable fit could be 
given. The width of the distribution is found to increase strongly 
with redshift. This is clearly a consequence of the increasing 
density inhomogeneities. Overall we find a large uncertainty of 
the proximity effect in individual sight lines. This is due to the 
limited w-scale coverage of the simulated spectr a, that are much 
smaller than in the observations analysed in Dall'Aglio et al.l 
d2008ah . 



9. Conclusions 

We have performed high resolution Monte-Carlo radiative trans- 
fer simulations of the line of sight proximity effect. Using dark 
matter only simulations and a semi-analytic model of the IGM, 
we applied our radiative transfer code in a post processing step. 
The aim of this paper was to identify radiative transfer influ- 
ence on the proximity effect for QSOs with different Lyman 
limit luminosities, different host environments, and at different 
redshifts. 

Due to the low optical depths in the Lya forest, we demon- 
strated that the on-the-spot approximation is insufficient in treat- 
ing the proximity effect. Diffuse radiation from recombining 
electrons significantly contribute to the ionising photon budget 
over large distances and need to be taken into account in the 
simulations. 

In the radiative transfer simulations we identified Lyman 
limit systems that cast shadows behind them in the overioni- 
sation profile. Any Lyman limit system between the QSO and 
the observer will result in a suppression of the proximity ef- 
fect. Behind such a system, the Lya forest is not affected by 
the additional QSO radiation and a detection of the proximity 
effect is not possible anymore. However, Lyman limit systems 
appear only on rare lines of sight and influence only marginally 
the proximity effect statistics. 

Mock spectra have been synthesized from the simulation re- 
sults, and similar methods as used in observations were applied 
to extract the observable proximity effect signature. We com- 
pared these results with the widely used assumption of the BDO 
model that only geometric dilution governs the over-ionisation 
region around the QSO. We have also used a random absorber 
model applied extensively in the analysis of observational data. 
Furthermore a semi-analytical model was employed which ana- 
lytically modelled the proximity effect signal onto spectra drawn 
from the UVB only simulations. 

We have shown that differences in the shape of the UVB 
spectral energy distribution function and the spectrum of the 



16 



Adrian M. Parti et al.: Cosmological radiative transfer for the LOS proximity effect 



QSO introduce a bias in the proximity effect strength. In the orig- 
inal formulation of the BDO model, the SED of the two quan- 
tities was assumed to be equal. Differences in the SED of the 
QSO and the UVB result in a weakening or strengthening of the 
proximity effect. This bias can be analytically modelled through 
introducing a correction term in the BDO formalism. We have 
tested the performance of the correction term. The correction 
term accounts mostly but not completely for the effect. However 
the remaining difference lies within the precision of the determi- 
nation of the strength parameter. 

In all simulations and models, we have confirmed a clear 
proximity effect profile. However, we found a strong scatter be- 
tween different lines of sight. Cosmological density inhomo- 
geneities introduced strong fluctuations in the normalised opti- 
cal depths of the proximity effect. They are caused by individual 
strong absorbers along the line of sight. The influence is signif- 
icant near the QSO and can dominate over the standard prox- 
imity effect profile if the QSO resides in a massive halo. The 
fluctuations in the normalised optical depth are smaller, if more 
luminous QSOs are considered. There the real space coverage of 
each log co bin is larger and therefore the fluctuations decrease. 
Additionally the fluctuations decrease with increasing redshift. 

We have studied the influence of different QSO host environ- 
ments on the proximity effect. The QSO was placed in the most 
massive halo, in a random filament, and in a random void. The 
local host environment of up to r ~ 2 Mpc h~ x is responsible for 
deviations from the BDO model near the QSO. However it does 
not affect the proximity effect at radii larger than r ^ 2 Mpc h~ l , 
and the smooth proximity effect signature was regained. 

In the three environments selected for this study, random 
density enhancement on scales up to r ~ 15 Mpc h were 
present. These overdensities introduced a bias in the mean anal- 
ysis, that decreases the proximity strength. It is a strong effect, 
seen in averaging the profiles over the 500 lines of sight, and then 
estimating the mean profile. This bias will affect the UVB mea- 
surement with this method. The influence of the bias decreases 
with increasing QSO luminosity. 

Our radiative transfer results are complementary to the anal- 
ysi s of the proximity effect zone with a semi-analytical model 
bv iFaucher-Giguere et al.l (|2008). There the authors find a sig- 
nificant dependence of the combined proximity effect on the 
QSO host halo both by the action of the absorber clustering 
and the associated gas inflow. Coinciding with our analysis, a 
stronger proximity effect signal is found at lower redshift due to 
the higher gas density inhomogeneities. The authors also discuss 
the nearly log-normal probability distribution of the flux decre- 
ment, and they find a strong dependence of the average overion- 
isation in the near zone around the QSO host from the halo mass 
on scales up to 1 Mpc proper radius. 

The proximity effect strength distribution derived along sin- 
gle lines of sight does not show such a dependence on the large 
scale overdensities. The distributions were always consistent 
with the random absorber model. This supports the analysis of 
real QSO spect ra and the resulting unbi ased estimate of the UV- 
background in lDalFAglio et al.l (l2008al) . We found that the dis- 
tributions are slightly shifted to negative log a values, since the 
co range covered in the simulation was smaller for the high lumi- 
nosity QSO. The fluctuations in the normalised optical depth on 
single lines of sight are responsible for a large scatter in the de- 
rived strength distribution. The scatter increases with decreasing 
redshift. Since the co scale in our analysis is limited, the fitted 
strength parameters show large uncertainties. This contributes 
to the width of the strength parameter distribution. A complete 



coverage of the proximity effect region will reduce these uncer- 
tainties on the strength distributions. 

To reliably study all the influences discussed here, the full 
proximity effect region needs to be simulated for which we have 
to make our radiative transfer code more efficient. It should also 
be noted, that Helium was not included in our simulations, but it 
could play an important role in additionally softening the QSO 
spectra. Helium would absorb a fraction of the hard flux pro- 
duced by the QSO and could increase the influence of the QSO 
SED. 

Finally we conclude that the QSO host environment, i.e. 
whether it sits in a halo, a filament or a void, influences the prox- 
imity strength only locally. Except shadowing by Lyman limit 
systems, the proximity effect does not show radiative transfer 
influence other then the SED bias. 

Appendix A: Strong absorber contamination of the 
normalised optical depth 

In Section|7]we have discussed the single line of sight variance of 
the normalised optical depth £ and its dependency on luminosity 
and redshift. In this Appendix we illustrate where and how this 
variance arises. 

The normalised optical depth £ is the ratio of the effective 
optical depth in the QSO spectra measured in log co bins over the 
effective optical depth present in the unaffected Lyo- forest. If the 
optical depth in a log co bin is larger than in the unaffected forest, 
£ will be larger than unity. The range covered by the log co bins 
in the spectra decreases towards the QSO. Therefore any strong 
absorption system present in the small ranges covered by high 
log co bins will bias the optical depth in the bins and therefore 
increase £ to values above unity. This bias gives rise to the large 
variance in £ at large log co. It is interesting to note, that an in- 
crease in the QSO Lyman luminosity by one dex also increases 
the covered proper distance in a log co bin by one dex. Thus the 
influence of strong absorbers is reduced, the more luminous the 
QSO is. 

In Fig. lA.ll we show five spectra taken from the z — 3 QSO 
situated in a void with a luminosity of L Vll = 10 31 erg Hz 1 s _1 . 
The spectra are plotted in log co space and the normalised optical 
depth measured in the log co bins is overplotted by solid dots. The 
spectrum at the bottom panel shows an uncontaminated prox- 
imity effect profile. However all the others panels illustrate the 
influence of high column density systems on the profile. In the 
second and third spectrum prominent absorption features elevate 
£ up to £ = 2. 

The distribution of f in a log co bin arising from the influ- 
ence of strong absorbers varies from bin to bin. In Fig. IA.2l we 
show the £ distribution in three bins for QSOs located in a void 
and a filament at redshift z = 3. Overplotted are the mean val- 
ues and the corresponding standard deviations shown in Fig. [12] 
We note that for the bright QSO the log co = 0.5 distribution is 
incomplete, since not all lines of sight reach the proper distance 
needed for this bin. 

The distributions are anisotropic with an extended tail to 
large optical depth. This behaviour reflects the probability dis- 
tribution of the simulated density field. The distribution of the 
high luminosity QSO is narrower than the one of the low lumi- 
nosity QSO. This is again due to the larger distances covered by 
the log co bins in the high luminous case. The width of the dis- 
tribution covers a large £ range. The contamination of a single 
line of sight by strong absorption systems is thus an issue that 
can only be eliminated by analysing large number of QSO sight 
lines to reliably determine the mean proximity effect profile. 
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Fig. A.2. Normalised optical depth £ distribution of the void (left panel) and filament QSO (right panel) at redshift z — 3 for three 
log a) bins. The upper panels show results for the L VhL = 10 31 erg Hz 1 s _1 QSOs, the lower one for the L Vll - 10 32 erg Hz 1 s _1 
QSO. The grey data points give the average normalised optical depths with their variance as shown in Fig.fTSl 
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Fig.A.l. Sample spectra taken from the z = 3 low luminosity 
QSO located in a void showing the radiative transfer simulations 
including the proximity effect (solid line) and without it (dashed 
line). The spectra are given in log u scale and for better visibility, 
no noise has been added to the spectra. The normalised optical 
depths £ at each log co bin is overplotted by solid dots to show 
the effect of absorption features on The dashed vertical line 
marks the start of the unreliability region of log co > 2.4. 
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